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Most calculations in cosmological perturbation theory, including those dealing with the infla- 
tionary generation of perturbations, their time evolution, and their observational consequences, 
decompose those perturbations into plane waves (Fourier modes). However, for some calculations, 
particularly those involving observations performed on a spherical sky, a decomposition into waves of 
fixed total angular momentum (TAM) may be more appropriate. Here we introduce TAM waves — 
solutions of fixed total angular momentum to the Helmholtz equation — for three-dimensional scalar, 
vector, and tensor fields. The vector TAM waves of given total angular momentum can be decom- 
posed further into a set of three basis functions of fixed orbital angular momentum (OAM), a set 
of fixed helicity, or a basis consisting of a longitudinal (L) and two transverse (E and B) TAM 
waves. The symmetric traceless rank-2 tensor TAM waves can be similarly decomposed into a basis 
of fixed OAM or fixed helicity, or a basis that consists of a longitudinal (L), two vector (VE and 
VB, of opposite parity), and two tensor (TE and TB, of opposite parity) waves. We show how 
all of the vector and tensor TAM waves can be obtained by applying derivative operators to scalar 
TAM waves. This operator approach then allows one to decompose a vector field into three covariant 
scalar fields for the L, E, and B components and symmetric-traceless-tensor fields into five covariant 
scalar fields for the L, VE, VB, TE, and TB components. We provide projections of the vector 
and tensor TAM waves onto vector and tensor spherical harmonics. We provide calculational detail 
to facilitate the assimilation of this formalism into cosmological calculations. As an example, we 
calculate the power spectra of the deflection angle for gravitational lensing by density perturbations 
and by gravitational waves. We comment on an alternative approach to CMB fluctuations based on 
TAM waves. An accompanying paper will work out three-point functions in terms of TAM waves 
and their relation to the usual Fourier-space bispectra. Our work may have applications elsewhere 
in field theory and in general relativity. 



Contents 

I. Introduction [U 

II. Notation d 

III. Scalar fields 

A. Plane waves 4 

B. Total-angular-momentum waves |4 

IV. Vector fields 

A. Plane waves 

B. TAM Waves 

1. The orbital-angular- momentum basis 

2. The longitudinal/transverse basis 

3. The longitudinal/transverse basis in terms of derivative operators 

4. The helicity basis 

5. Projection onto Vector Spherical Harmonics 

6. The plane wave expansion for vector fields 

7. Expansion of vector fields and power spectra 

V. Symmetric Tensor Fields 

A. Introduction and Plane Waves 

B. TAM waves 

1. The orbital-angular- momentum basis 

2. The longitudinal/vector/transverse-traceless basis 

3. The longitudinal/vector/transverse-traceless basis in terms of derivative operators 

4. Summary and helicity basis 



2 



5. Projection onto Tensor Spherical Harmonics IS 

6. Plane-wave expansion for traceless tensor fields 21 

7. Expansion of tensor fields and power spectra [22 

VI. Calculation of Lensing Power Spectra 23 

A. Scalar metric perturbation 24 

B. Tensor Metric Perturbations [24 

VII. Boltzmann equations for CMB Fluctuations 
VIII. CONCLUSIONS 

Acknowledgments [28| 

A. Divergence of the vector harmonics 28 

1. Gradient of scalar TAM waves 28 

2. Divergence of vector TAM waves [29 

B. Divergence of tensor harmonics [3(3 

C. Irreducible tensors [3f| 
References 



I. INTRODUCTION 



Much of modern cosmology involves the study of the origin and evolution of scalar, vector, and tensor fields. 
Examples of scalar fields include the inflaton [l| and the quintessence field @. Perturbations in the inflaton are 
considered as seeds for primordial perturbations to the curvature, and there is active investigation of the effects 
of quintessence perturbations on the evolution of density perturbations at late times. Magnetic fields provide an 
example of vector fields in cosmology 0, 0] , and vector fields have appeared elsewhere as well Q . The most general 
perturbation to the spacetime metric involves a rank-2 tensor field [6( , the six components of which, as is well known, 
can be decomposed into a trace, a longitudinal mode, two vector modes, and two transverse-traceless modes, the latter 
of which propagate in general relativity as gravitational waves. A stochastic background of cosmological gravitational 
waves may have been produced during inflation Q and are now being actively sought through the curl component (B 
mode) they induce in the anisotropy of the cosmic microwave background (CMB) polarization 

The vast majority of the literature on such fields and perturbations proceeds by decomposing the perturbations 
into Fourier modes (or plane waves) e lk x , each of which then evolves independently to first order in perturbation 
theory A primordial random field, such as that produced by inflation, is then assembled by adding all such plane 
waves with the Fourier amplitude for each wavevector k selected from a Gaussian distribution with a variance given 
by the power spectrum P(k). 

However, observations of the Universe are performed on a spherical sky. It may thus be advantageous, in some cases, 
to consider decomposition of the scalar/vector/tensor fields under consideration into a basis that reflects the rotational 
symmetry of the spherical sky. With this motivation in mind, we introduce here total-angular-momentum (TAM) 
waves for scalar, vector, and tensor fields. We provide a complete orthonormal set of basis functions for scalar, vector, 
and tensor fields on three-dimensional Euclidean space, of fixed orbital angular momentum. These basis functions are 
cigenfunctions of the Laplacian operator. There are three vector TAM waves for each total angular momentum, and 
we decompose these three into a basis of fixed orbital angular momentum (OAM), a basis of fixed hclicity, and a basis 
that separates the longitudinal (L) and two transverse modes, E and B, of opposite parity. There are similarly five 
sets of TAM basis functions for traceless rank-2 symmetric tensors, and we provide similar decompositions into three 
sets of bases: an OAM basis, a helicity basis, and a basis that decomposes into a longitudinal (L) TAM wave, two 
vector waves, VE and VB, of opposite parity, and two transverse waves, TE and TB, of opposite parity. A random 
field can then be assembled by adding all such TAM waves, with the amplitude for each TAM wave of wavenumbcr 
k selected from a Gaussian distribution with variance P(k), as will be detailed below. 

The L/E/B basis for vectors, and the L/VE/VB/TE/TB basis for traceless symmetric rank-2 tensors, are first 
derived from the OAM basis. The helicity bases are then simply related to these L/E/B and L/VE/VB/TE/TB 
bases. We then present an alternative derivation of these TAM-wave bases by introducing sets of vector and traceless- 
tensor differential operators that, when applied to the scalar TAM waves, yield the L/E/B and L/VE/VB/TE/TB 
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TAM waves. This operator approach then allows one to write an arbitrary vector field V a (x) in terms of three 
covariant scalar functions V L (x), V E (x), and V B (x) and an arbitrary tensor field /i a fc( x ) m terms of five covariant 
scalar functions /i L (x), h VE (x), h VB (x), h TE (x), h TB (x). The operator approach also allows one to obtain sets of 
vector and tensor basis functions from any other set of scalar basis functions. 

We also provide the projections of the TAM vector and tensor waves onto vector and tensor spherical harmonics. 
This is equivalent, as will be seen below, to providing the n, 9, and (j> components of the TAM vector and tensor 
waves. The utility of TAM waves, as well as these projections, is illustrated with a re-derivation of the power spectra 
for weak lensing by density perturbations and gravitational waves. 

The bases we provide here for three-dimensional fields are to be contrasted with earlier work 0] , developed for CMB 
polarization and weak- lensing shear, on bases for tensor fields on the two-sphere S 2 , and with bases for vector fields 
(the weak- lensing deflection angle) Q ; we will show below, though, how the 9-6 components of the three-dimensional 
TAM waves map onto the familiar vector/tensor spherical harmonics. Ref. [10Ml2| present the E/B modes of three- 
dimensional vector and tensor harmonics in open and closed Friedmann-Robertson- Walker space. The TAM-wave 
basis for scalar fields has already been employed in cosmology [l3l - [l7l | (sometimes referred to as a "spherical- wave" 
or "Fourier-Bessel" expansion) . The vector TAM waves are familiar from electromagnetic theory [l9|, H3] . Some steps 



along these lines for tensor fields were taken in Ref. [21[ , although they retained plane waves for the spatial dependence. 
There are some resemblances to Ref. [22j| . who were considering classical cosmological tests. Rcfs. [23|, [24| have taken 
significant steps in the direction we pursue here for the description of weak lensing by density perturbations, and 
there are some analogues to this work in the gravitational- wave literature (see, e.g., Ref. j25|). 

Below we begin in Section |TT] with a brief discussion of our notation. Section UTTl presents TAM waves for scalar fields 
beginning, by way of introduction, with a review in Section UlIAl of the Fourier expansion of scalar fields. Section ITVl 
discusses vector fields, beginning in IIV Al with plane waves and moving on in HVBI to TAM waves with vector fields. 
There the TAM waves of fixed total angular momentum are decomposed into OAM, L/E/B, and helicity bases. We 
introduce in Section TlV B 31 a set of derivative operators that, when applied to the scalar TAM waves, provide TAM 
vector waves in the L/E/B basis. We also show here how this operator approach can be used to find scalar functions 
associated with the L, E, and B components. The rest of Section HVl discusses the projection of the TAM vector 
waves onto vector spherical harmonics ([IV B 5p . results that are useful, e.g., for observational quantities like the lensing 
deflection field that are represented as vectors on the two-sphere; the transformation between vector plane waves and 
vector TAM waves QIVB 6|) ; and the expansion of vector fields in terms of TAM waves and the relation between the 
TAM-wave power spectra and the more familiar plane- wave power spectra (jIVB 7[) . Section [V] provides a discussion 
of tensor TAM waves with an organization that parallels precisely that for vector waves in Section IIVI Section IVII 
presents, as an example of the utility of the TAM-wave formalism, a calculation of the power spectra for the deflection 
angle from gravitational lensing by density (scalar) perturbations and gravitational waves (transverse-traceless tensor 
perturbations) . Section IVIII discusses the prospects for writing the Boltzmann equations for the evolution of CMB 
fluctuations using the TAM-wave formalism. Section IVIIII provides closing remarks. Appendixes [A] and |B] provide 
calculational details, and Appendix [C] provides a proof that the functions obtained by the action of irreducible-tensor 
operators on TAM waves are TAM waves of the same total-angular-momentum quantum numbers JM. 



II. NOTATION 



Throughout this paper we use the symbol ^"(x) to denote solutions to the Hclmholtz equation, 

(V 2 +fc 2 )* fc (x) =0, (1) 

where x is the spatial position, and k is the magnitude of the wavevector. In order to reduce clutter, we will often 
suppress the k superscript. We will be dealing with solutions ^^ JM ) (x) that are eigenstates of total angular momen- 
tum and its z component labeled by eigenvalues J and M, respectively. We will also obtain scalar, vector, and tensor 
solutions to the Helmholtz equation, and we will denote those solutions (actually, the components of those solutions) 
by ^(ja/)(x), $(jm)oW, and v I'(jA/)a6( x ) (where we have suppressed the k superscript, as we will do frequently 
throughout), respectively. The number of indices in the subscript, outside the parentheses, indicates whether the 
quantity is a scalar, vector, or tensor. As we will see, the vector and tensor eigenfunctions of fixed JM can be decom- 
posed into states of fixed orbital angular momentum, fixed helicity, or a longitudinal/transverse decomposition. These 
will be labeled by a superscript. For example, the vector eigenstate of total angular momentum JM for wavevector k 
with orbital angular momentum I will be ^j'j^^x), and the vector TAM waves in the transverse/longitudinal basis 

will be referred to by *&(jM)a( x )i ^ or a = L,E,B, and in the helicity basis by \I/^^ a (x), for A = 0,±1. Again, the 
k superscript will often be suppressed. We often refer to V a as a "vector," although strictly speaking, it is the dual 
vector associated with the vector V a ; there should never be any confusion given that the dual vector has a lowered 



4 



index and the vector a raised index. The indices are raised and lowered with a metric g a b, and the antisymmetric 
tensor is e a b c - Since we are dealing with flat three-dimensional space, the metric may be taken to be a Kronecker 
delta with Cartesian coordinates, in which case the raising and lowering of indices is trivial. However, we will at times 
work in spherical coordinates r, 6, <fi in which case g a b is not trivial. In some places we will deal with functions on the 
two-sphere § 2 , and in these cases we denote the metric and antisymmetric tensor for the two-sphere by qab and £ab , 
respectively, with capital indices. 

We will also, by way of introduction, deal with plane-wave solutions to the Helmholtz equation. We will label the 
scalar, vector, and tensor solutions by ^ k (x), ^(x), and ^^(x), respectively. An additional superscript will denote 
the decomposition into OAM, helicity, or longitudinal/transverse eigenstates. 1 For reference, we list in Table I the 
symbols used in this paper. 



III. SCALAR FIELDS 



A. Plane waves 



We begin with scalar fields to provide a simple introduction. Our aim is to find solutions 0(x) to the scalar Helmholtz 
equation, (V 2 + fc 2 )0(x) = 0. The most general solution can be written in terms of plane waves 4' k (x) = e jk x , 
cigenfunctions of the momentum operator — iV. The set of solutions for all k constitute a complete orthonormal basis 
for scalar functions 0(x), normalized so that 

d 3 x* k (x) [* k '(x)j* = (27r) 3 (5 D (k-k'), (2) 

where <5u(k — k') is a Dirac delta function. The most general scalar function can then be expanded, 

k 0(k)* k (x), where «£(k) = f d 3 x0(x) [* k (x)] * . (3) 



(2tt) 3 

The power spectrum P(k) for a scalar field is then defined by 

;<Kk)0*(k')) = WMk- *)p(k), (4) 

where the angle brackets denote an expectation value over all realizations of the random field. 



B. Total-angular- momentum waves 



Our aim here, though, is to find solutions that are eigenfunctions of angular momentum. This is easily done with 
the plane- wave expansion, 

e ik ' x = ^4m l ji(kr)Y$ m) $)Y Qm) (h), (5) 
Im 

where ji(x) is a spherical Bessel function, and Yn m )(fi) are (scalar) spherical harmonics. 2 We then find that if we 
choose total-angular-momentum (TAM) basis functions, 

*f M (x)=i i (fcr)F (im) (fi), (6) 
where r = |x| and n = x/r, then an arbitrary scalar function can be expanded as 

= Hi J^^(lrn){k)Am l ^\ lm) {^), (7) 

Im 



1 Note that there is no such superscript for scalar waves, as the OAM and TAM waves coincide for scalar fields, since they have no spin. 

2 We choose the spherical Bessel function of the first kind ji(kr) rather than the second kind n;(fcr), so that the TAM waves are regular 
at the origin. There may be cases, for example in application of this formalism to emission or scattering of gravitational radiation, in 
which the second function n;(fcr) may need to be introduced. 
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x = x, y,z, r = |x|, and n 


a point in , its norm, and a unit vector in its direction 




k, k a , and k — k 


Fourier wavevector, its components, and its magnitude 




k, e, 


etc. and k a , 6 a , etc. 


unit vectors in the k, 6, etc. directions and their components 




v, v„ 


covariant derivative wrt x and its components 




* and f 


as superscripts represent complex conjugation and hcrmitian conjugate 




(X) 


average over all realizations of random variable X 


(|4| 


S D (k-k') 


the one-dimensional Dirac delta function 


17TT 


fc(k-k') 


the three-dimensional Dirac delta function 


TTTT 


a, b, c, . . . 


three-dimensional tensor indices 




A,B,C,... 


two-dimensional tensor indices 




£« 


polarization vector 




5ij 


Kronecker delta 




g a b and e abc 


metric and antisymmetric tensor in IR 




gAB and eab 


metric an antisymmetric tensor on a 


1)891) 


e ab 


polarization tensors 






polarization tensor for tensor plane wave 


(I96|) 


h a b{k) 


Fourier transform of /i a &(x) 


m 


h s (k) 


amplitudes for tensor plane-wave components 


(|67p 


L and L a 


orbital-angular-momentum operator and its components 




S 


spin operator 




J and J a 


total-angular-momontum operator and its components 




J and M 


Quantum numbers for total angular momentum and its z component 




I and m 


Quantum numbers lor orbital angular momentum and its z component 




1 


' Ii h h \ 
y mi m 2 m,3 J 


Wigner-3j symbol 




1 


h h h 1 

mi 7712 U13 J 


Wigner-6j symbol 


(iBll 


(Zimi/ 2 m 2 JM) 


Clebsch-Gordan coefficient 




(1201) 


V a (x) 


vector field 


CUT 


ftai>(x) 


tensor held 


H67J) 


y Q (k) 


Fourier coefficients for vector field 


GUI 


V™(x) 


scalar fields for a = L,E, B components of vector field 


(|64|) 


^(x) 


scalar fields for a = L, VE, VB, TE, TB components of traceless tensor field 


((TW 


P(fc) 


power spectrum for density perturbations 


o 


P L (k),P T (k) 


i c i • j l ■ l i i l c ir"ii 

power spectra tor longitudinal and transverse modes ol vector field 


dizj) 


P+(k),P-(k) 


power spectra for left- and right-circularly polarized vector fields 




P h (k) 


power spectrum for transverse-traceless mode of tensor field 


(loci) 


* fc (x) 


Solutions to the Helmholtz equation for wavenumber k 


/TTK 


0(x) 


scalar functions 




0(k) 


the courier transform ol </>(x) 


© 




1 AM- wave transform ol <p(xj 




f(x) 


scalar for longitudinal component of h a b 




to a (x) 


transverse- vector field for vector component of hab 






transverse-traceless part of h a b 


HSU) 


* k (x) 


scalar plane-wave mode 




*2' k (x) 


i i l r l ■ j ■ t i c\ 

vector plane-wave mode lor polarization a = L, 1, 2 


in 


*£'*(x) 


circularly-polarized vector plane- wave mode lor helicity A = ±1 




j; fx) and n ; (a;) 


spherical Bessel functions of the first and second kind 






scalar spherical harmonic 


doll 


e m 


Qnhpripfil bfiQit: fnv vpptnr 




i m 


spherical basis for tensor 






vector spherical harmonic of OAM for / = J — 1, J, J -\- 1 


((20)) 




vector spherical harmonic in the longitudinal/transverse basis for a = L, E, B 

i O / i i 


(|44p 


*rW( A ) 


vector spherical harmonic of helicity A = 0, ±1 


CD 


^ / (ji\/)t t i)(") 


tensor spherical harmonic of OAM for I = J — 2, . . . , J + 2 


CD 




tensor spherical harmonic for a = L, VE, VB, TE, TB 


([92]" 




tensor spherical harmonic of helicity A = 0, ±1, ±2 


([99)1 



(1771 ) V / 


scalar TAM wave 




^ J Ai J a V / 


vector TAM wave of OAM Z, for / = J - 1, J, J + 1 


(I2U1) 

1 V 


*7™n fx) 

( J Ai )a \ J 


vector TAM wave in the longitudinal/transverse basis for a = L,E,B 


J35T 


fx) 

( J Ai ) a V / 


vector TAM wave of hclicity A = 0, ±1 


(1361) 


t// ya yA 

( J Ai ) ' ( J Ai ) ' (J 111 ) 


vector TAM-wave coefficients 


(El) 


*5 i.n i,(x) 

( J Ai ) a o \ / 


tensor TAM wave of OAM Z, for I = J - 2, . . . , J + 2 


(1711) 




tensor TAM wave for a = L, FS, V\B, T_E, T£ 


(JUT 


^ (JAf)a(A^V 


tensor TAM wave of helicity A = 0, ±1, ±2 


dHTb 


n (JM)> "(JM)l "(JA/) 


fpnQnr TATYT-wrtvp OYntmsinTi ppipffipipntt; 

L dloWl J_ A1V1 VV CL V 1_jjY UclllolWll LjvJ1_^11H_.11_i11Lij 


(|102|) 


■^fJM)) rt (JM)l ^(JAf)' 


pil n p-\n7Fi^7"p pYnan qi rin priPTn pi p~n1~Q tciv vppt pi i~q 

UldllC VVaVC CA |Jdiloll_Hl L^UClllL^lCllLO 1U1 V CL LU1 O 


(LiUI 


p/ pa pA 

a (JM); a (JM)> D (JM)> 


plane-wave expansion coefficients for tensors 


(IUDII 


JJ a , M a , ana K a 


differential operators 


/lopiii 
(loUII 


TV A/f , nnrl £T 


nnQT 1 o "T" Ave n Arm o 1 ( 1\J \ f oiirrrMTi"! q 1 ^ 7\ ■■'/" . /V ^ + a 
(jpt^IclLLHo I1UI Illtii [1\ a J LajIl^cllLlcll wV'iJ^ a , 1 LU 


/Roll 




i fl"^T"i^Ti "f" i o 1 *" I / i >' o "f - /~i»""c rrori Drof in rr \Tf ^ "Fv/^m \T/ ^ T ~ r 




W a 


differential operators generating ^' / ( < jM)a£> from 




'XI J rri .J rri J 

J-aU j Aq) J-Xl 


transtormation matrices between dinerent vector bases 


13JJ 


U a i, Ux a , Uxi 


transtormation matrices between different tensor bases 


((88)) 


A Q 


lensing deflection field 


dii7| 


n ab 


projection tensor 


(I38|) 


and 770 


conformal time and its value today 


(|117|) 


<£>(n) and S2(n) 


projected lensing potentials 


(|118|) 


f(.IM) anu H(JM) 


1 pn qi n o'-Pinf PYifi a 1 qpi ViPTi pa 1- h a rm r>n i p poptti pi pti^q 

ICllOlll cL \J \J LrOllLjlCLl BU11C1 1 >_ til 11 Cl.I IIUJIIIL IjULlllLlLlllo 


u ^ ^ f u 


C vv and CV si 


1ptigJt-ip" fincnilar nnwPT m^ppfrfi 

ICllolllt Ctllti Ll±ClJ_ yj \J VV V^L DL/Cl^ul Oi 


(|119|) 


T sca (fc) 
v / 


density-perturbation transfer function 


(11211) 


(2(7^) and D\{rf} 


scale factor and density-perturbation growth factor 


(|12ip 




primordial amplitude for TAM wave of density perturbation 


(|121f 


' (A) 


lensing transfer function for density perturbations and gravitational waves 


(11231) 


Uj ana Oj 


lensing angular power spectra 


ii or 

(|12o|) 




^~ittip fii rpi hif i nn tpit rri - Ci ~\n i~ Ci ~t~ i r^n si 1 ursvp 

LIIIHj \_< V L*l LI L 1W11 til d V 1 Lid/LilWllCll WdVC 






roHiQ 1 di rv/^n tii npf i An c tat" n*T" o^nf oliAn o 1 AiroirA f ran offir fun i An o 
IdtlleLl cl^cillLlIlLLllJIlb 1(JI ^1 dvlLdLlUIlcll- WtiVt; Lldllblcl ILlIlL-LKJIlh 


HI ^9h 


"V*-! q, //j 


ran i ni~ i nil PiPT*ti i rnatinn 

1 dUld L1LJ11 l/Cl L HI U (XL 1U11 


(11 4-311 




expansion coefficients for radiation perturbation 


111 /I 


•=k,JM, a n 

-n' ( x >q) 


1AJV1 eigentunctions of x and q 


(1143 P 


q 


direction of photon momentum 


11 1 A 9 h 
(|14d|) 




mth component of irreducible-tensor operator of rank I 


|C2| 




vector operator for a — L, E, B 


<ETJ 


X>' , 

mm' 


Wigner rotation matrix 




TZ 


an 0(3) rotation 


(IC2l) 



TABLE I: A list of mathematical symbols used. The number in the right-most column indicates the equation in or near 
the symbol is first defined or used. 



(ft 



W*f im) (x) 0(x)= / d 2 k0(k)r ( t m) (k). 



with 

4>(im)(k) 

Here we have used the relations, 

J k 2 dkji(kr)j l (kr>) = ^6 D (r-r'), £ F (te) (n)^^ (n') = S D (n - n') 

The orthonormality relation for the basis functions is 

(27T) 3 



where 5« is the Kronecker delta. The basis functions also satisfy, 



-5 D {k - k'), 



E 



p dk 

(2^)3 



47ri'*f Zm) (x) 4^*f im) (x') =£n(x-x') 
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which demonstrates that the ^% m ) (x) constitute a complete basis for scalar functions on M 3 . The products of TAM- 
wave coefficients have expectation values, 

(^i m) {k)4>tv m ')(k')) = ^rS D (k - k')8 w S mm >P(k). (12) 



IV. VECTOR FIELDS 
A. Plane waves 

We now generalize to vector fields. Again we begin by reviewing plane-wave vector solutions to the Helmholtz 
equation. Three solutions to the vector Helmholtz equation, (V 2 + fc 2 )4' a (x) = can be obtained, for each Fourier 
wavevector k, as 



v^< k (x) = (lA-)V„$ k (x)=ie 



^• k (x) = -^V h ^ k (x) = — ^(V a V b -. 9ah V 2 )£ b * k (x) = iP^L^ke*-, (13) 
k «|k x z| fc|k x z| 

where z is a unit vector in the z direction, and e a b c is the totally antisymmetric tensor. Here v I'^' k (x) is a longitudinal 
vector field, and Vf'*' k (x) and '^'''(x) are the two linear polarizations for the transverse part of the vector field. We 
could have written Eq. ([13"]) more simply as ^^' k (x) = ie„ (k)e x , with = k a and e\' 2 two other unit vectors 
orthogonal to k a and to each other. We have written in Eq. (|13[) one choice for these polarization vectors explicitly 
in terms of a fixed unit vector z to motivate a choice of polarization vectors for the TAM waves later. 
These mode functions are normalized so that they constitute a complete orthonormal set, 



J d 3 ztf"> ka (x) [** k (x)]* = (27r) 3 Mk-k')<W, 



(14) 



where a, /3 = {L, 1,2}. The three mode functions are, furthermore, orthogonal at each point. An arbitrary vector 
field V a (ii) can then be expanded as 

K(x) = / [^(k)^< k (x) + y J (k)^ k ( x ) + ^ 2 (k)* 2 - k (x)] , (is) 

in terms of Fourier expansion coefficients, 

V L {k) = J d 3 xV a ( X ) [^ k (x)]* = -J d 3 x [* k (x)]* iv a K(x), 

V*(k) = J A-ffx) [^' k (x)]* = J d 3 x [*k( x) ]*_l_ eahc ^vV c (x), 

V 2 {k) = I d 3 x V a (x) [* 2 < k (x)] * = J d 3 x [* k (x)] * ( Va V 6 - 5afc V 2 )F b (x). (16) 

We obtain the last equality in each line by integrating by parts. 

The statistics of the vector field are given in terms of power spectra Phik) and Pr{k) for the longitudinal and 
transverse components, respectively, that satisfy 



V L {k)V L *(k')) = (2ir) 3 6 D (k-k')P L (k), 
^CkjV^k')) = {2*) 3 5 D {k-k')P T {k), 

V 2 (k)f 2 *(k')\ = (27r) 3 5 D (k-k')P T (k). (17) 



The decomposition of the transverse component into the two modes V a (a = 1, 2) is not rotationally invariant — the 
decomposition would be different if we had chosen a different direction for z — so the power spectra for the two must 
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be the same. However, we can alternatively decompose the transverse- vector modes into plane waves of right and left 
circular polarization, or positive and negative helicity, 

vl/±' k (x) = -i=(^ k (x)±*vl/^ (x) ). (18) 

Since ^^' k (x) = — (i/fc)e a b c V fa 4' 1,k,c (x), these modes are invariant under rotations about the k direction and thus in 
some sense more "physical" than the 1 and 2 linear polarizations. It is possible (although it would require parity 
breaking) that P+(k) and P-(k) could differ. In the absence of parity breaking P+{k) = P-(k) = Pr(k). 



B. TAM Waves 



The aim now is to find vector- valued functions V a (x) that satisfy the vector Hclmholtz equation, (V 2 + fc 2 ) V a (x) = 
0, for definite wavenumber magnitude k, and that transform under spatial rotation as representations of order J. In 
other words, we seek eigenfunctions of total angular momentum J = L + S, where L a = —ie a b c x b V c is the orbital 
angular momentum and S is the 5 = 1 spin associated with the vector space spanned by a set of basis vectors at each 
spatial point. This differs from the case of scalar fields where, with spin 5 = 0, total-angular-momentum eigenstates 
coincide with orbital-angular-momentum states. 

Our strategy will be to first construct vector- valued eigenfunctions of total angular momentum that are also eigen- 
functions of orbital angular momentum We will then construct linear combinations of states of definite total angular 
momentum JM that are curl- free (the longitudinal component) and divergence- free (the two transverse components). 
We will then decompose the TAM waves also into helicity eigenstates. 



1. The orbital- angular-momentum basis 

From the usual set of Cartesian basis unit vectors = 8%, for [i = x,y, z, one can construct a spherical basis e™ 
for m = +1, 0, —1, through 3 

,o _ z ± 



< = e' a , e±=T{el±zel)^2. (19) 

These constitute a complex but global basis, so these unit vectors commute with differential operators. Under spatial 
rotations, they transform as an I = 1 representation. We know for the spatial part of the eigenfunction that the 
conventional scalar- valued spherical harmonics Y(; m ) (n) form a representation of order / of the spatial rotation group. 
Vector- valued eigenfunctions of total angular momentum arc therefore constructed via the usual scheme for adding 
two angular momenta (20| . 

*(*iif)o( x ) = 3l{kr)Y l {JM)a {h) ee 5^<lmim|JM)i,(*r)y (lm) (fi)e*, (20) 

mm 

where (Zimi^^l JM) are Clebsch-Gordan coefficients. Here, J — 0,1,2,..., M = — J, — J + 1,..., J — 1, J, and 
I = J — 1, J, J + 1. The TAM waves ^'^^(x) are also eigenfunctions of orbital angular momentum squared 
L 2 = L a L a with eigenvalue 1(1 + 1). The angular parts Y^ JM ^ a (ii) (three-dimensional vector spherical harmonics of 
given total and orbital angular momentum) are normalized to 



d 2 n 



Y l {J ? MI) (h) Y l {JM)a {rL) = 8 lv 8jj,8 MM >. (21) 



There are three eigenfunctions for given total angular momentum JM distinguished by their orbital angular momen- 
tum /. The TAM waves are normalized so that 



d 3 x 



^tf; M) jx)] * 4 7 r l J '*;' J ^ ) (x) = 8u>5jj.6' MM ^8 D (k k'). (22) 



3 Barred indices like m are reserved for the order-1 spherical basis. 
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We also have that 



E 

JMl 



k 2 dk 



W^(x') 4™%* (x) =fe(x-x'), 



(23) 



which demonstrates that the ^(j^f) a ( x ) constitute a complete basis for vector functions on M 3 . To show this, we use 
the definition in Eq. (f20| to rewrite the left-hand side as 



k 2 dk 

JMl 

The k integral is 



k 2 dkji (kr)ji (kr 1 ) 



E fecfiO 



JMl 



Y(JM)a^)- 



- f k 2 dkj l (kr)ji(kr') = \s D (r-r') ) 



(24) 
(25) 



and the sum then becomes 



E[ y (W fi/ ) 



Y(jM)a^) = EEE (lmH^)(lmW|JM)[y (im) (n')e* a ]*y (w) (n)ef 

7Mi mm m'm' 

= EE E (l^m|JM)<lmW|JM)[y (lm) (fi')]*l r (Im0 (fi)* fiMft< 

JM! mm m'm' 

= E E [^Om) ("')]* y(/m')(n) (5 mr7 i "5mm' 
Imfh m'fh' 

= [Y { i m) (h')YY (lm/) (h) = 6 D (A-n'), 



(26) 



from which Eq. ([23)) follows. Note that this also demonstrates that the y/ JM -,(n) constitute a complete basis for 
three-dimensional vectors on the two-sphere. 



2. The longitudinal /transverse 

The next step will be to construct linear 
To do so, we must calculate the divergence a 

m is 



istruct linear combinations of these OAM waves that are longitudinal and transverse, 
j the divergence and curl of ^'^^(x). The result for the divergence, detailed in Appendix 

(-^f^(kr)jj(kr)Y {JM) (h), I = J - 1, 

rV a *[ JM)a ( x ) = I 0, I = J, (27) 

I - ^ jj( kr ) Y (JM) (&), Z = J + 1. 

Since the parity of a given OAM state is (—1)', and the basis vectors in Eq. (|20[) of odd 
vector fields of parity (— 1) J and (— 1) J+1 to be, respectively, 4 

yfjM)a(*)=MijM)a(x)> and *f«),W = 

The basis 



I parity, we choose 



i functions for the longitudinal field may then be taken 

.-r.T, / \ . J .-r.J-l I \ 



2J + 1 (JM)oW 



2J + l*( / " / ) a(x) 



transverse- 
• (28) 



to be 



*fjM)a( X ) = * 



(29) 



4 The parity of (— 1) J for the vector i? mode and (— 1) for the vector E mode differ by —1 from the parities of the E/B tensor spherical 
harmonics. However, the expansion coefficients for the vector E and B modes have, as we will see below, parities (— 1) J + 1 and (— l) -7 , 
as do the tensor-spherical-harmonic expansion coefficients. The reason traces back to the transformation property of the vector field 
under a parity inversion. 
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The prefactors have been chosen so that the three sets of eigcnfunctions are normalized as in Eq. ([22]) . We thus have 
a complete orthonormal set of basis functions, of fixed total angular momentum, for the transverse and longitudinal 
components of a vector field. 



3. The longitudinal/transverse basis in terms of derivative operators 



There is, however, an alternative and useful route to these longitudinal and transverse basis functions. In Ap- 
pcndix[C]it is proved that if an operator O is an irreducible tensor under rotations, then J 2 OY^ JM ) = J{J + l)OY( JM \ j 



and JyOYi 



(JM) 



operators to scalar TAM waves. 
Consider three vector operators 



MOY^jMy We can therefore construct vector TAM waves by applying appropriately defined vector 



D a = -V„, 

fc 



K a = -iL a , M a = e abc D b K c 



(30) 



These are irreducible- vector operators, and they all commute with V 2 . They therefore yield, when acting on scalar 
TAM waves, TAM vector waves of total angular momentum JM that are also solutions of the vector Helmholtz 
equation. These three sets of vector fields must be linear combinations of ^/ JM Jx), for I = J — 1, J, J + 1 or 

a.k 



j M j (x) for a = L,E, B. Since the three operators satisfy 

D a K a = K a D a = 0, D a M a = 0, M a D a 



2, K a M a = M a K a = 0, 



(31) 



it follows that D a generates the longitudinal vector field Da$>^ JM ^ (x) oc >I^J M ^ a (x), while K a and M a generate 
divergence- free vector fields. Since K a is axial- vector-like and M a vector- like, parity considerations tell us that K a , 
generates the B mode, K a ^ JM ^ (x) cx ^fr M ) a (x), while M a generates the E mode, M a ^^ JM ^(x.) oc \I/^^ a (x). 

The operators D a , K a , and M a are also operators in the Hilbert space of vector- valued fields, so we can calculate 
their hermitian conjugates to be 

(M ) t = ~M a + 2D a . 



Thus, when acting on the three operators have norms 

(D a yD a = 1, (K a )^ K a = L a L a = J(J + 1), (M a )Hl a = L a L a = J (J + 1). 



(32) 



(33) 



These results enable us to normalize the vector TAM waves and to reproduce the longitudinal/transverse basis. This 
operator approach has the advantage that many calculations involving vector or higher-spin TAM waves can be 
reduced to the algebra of operators that act on scalar spherical waves. The following properties of the three operators 
will be useful in calculations: 



[D ai D b ] = 0, [K a , D b ] = e abc D c , [M a , D b ] = g ab - D a D b , 
[K a , K b ] = e abc K c , [K a , M b ] = e abc M c , [M a ,M b ] = -e abc K c 



(34) 



We can gain insight into the operators D a , K a , and M a from the far- field limit kr — > oo, where D a is approximated 
by an ordinary vector in the radial direction, and K a and M a asymptote to two orthogonal vectors in the plane 
perpendicular to the radial direction, when they act on a scalar TAM wave. The factor of i in the definition of D a is 
chosen so that D a = — k a in this limit. The sign convention for the E/B vector TAM waves is chosen so that if we 
rotate the E mode by +90° about the direction of wave propagation we obtain a B mode. 
To summarize, the decomposition into longitudinal and transverse modes is 



B mode : * 



(JM) 



„(*) = 



E mode : ^ 



(JM) 



„(*) = 



(JM)\ 



(x) = ^ £Qbc V h *f JM) (x) = zvE-( 7 7M)a (x), 



(JM)( X ) — J. e abc 



V^f /A c f) (x) 



2J+ 1 



1/2 



J 



2J+1 



1/2 



J+l 



longitudinal mode : 



*fjM)«( X ) 



£>a*(JM)(x) = i 



J 



2J+ 1 



1/2 



JJ\/)a( X ) 

V2J+ 1 



1/2 



(35) 
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4- The helicity basis 



We can define another basis, denoted by the helicity A = 0, ±1, by 



1 



*fA/)» = ^ [*(W(x) ± i*f JM)o (x) 



(36) 



These are eigenstates of the helicity operator H = S • p, where (S b ) ac — i^abc is the spin operator and p a = — iV a /fc 
the normalized momentum operator, with eigenvalues A. 

We may summarize the transformation between the three bases — labeled by I = J, J— 1, J+l for the orbital-angular- 
momentum basis, a = B,E,L for the longitudinal/ trans verse basis, and A = 1,0, —1 for the helicity basis — by the 
transformation matrices, 



/I 

















J+l 

2,7+1 
J 

2,7+1 



2.7+1 
27+1 / 



— — 

tL = I o oi 

— — 

V2 V2 U 



v'2 





.7+1 
2(2,7+1) 

J 



7+TT 1 



2,7+1 
.7+1 



2(27+1) 
.7+1 



2J+1 



\ V2 V 2(27+1) V 2(2.7+1) / 



(37) 



5. Projection onto Vector Spherical Harmonics 

Here we have constructed three different bases for three-dimensional vectors on R 3 . We now show how the angular 
components project onto the more familiar vector spherical harmonics Y|j M ^ a (n) and Y^ M ^ a (h), for two-dimensional 
vectors that live on the two-sphere. These vector spherical harmonics are given by 

Y(jm)°( & ) = 77p^ry v±aY(JM)(fi) ' F (W( fi ) = ^j(7+i) ea6c " 6vcy(JM)(fi) ' (38) 

both of which have h a Y^j M - )a (h) = = n a Y^ M ^ a (ri). Here W± a = n a 6 Vb is the gradient operator in the 6-<j) 
space, and II a b = g ab — h a h b projects onto that space. In addition, we can define a third vector spherical harmonic 
y^Mjotn) = —n a Y(jM){^) to account for the component of a three-dimensional vector in the normal direction. 

This set of three vector spherical harmonics provides a complete set of orthonormal basis functions for three- 
dimensional vectors that live on the two-sphere. We can obtain these vector spherical harmonics using an operator 
approach that parallels that we developed for TAM waves. Define three dimensionlcss irreducible- vector operators, 

N a = -n a , K a = -iL a , M± a = e abc N b K c . (39) 

These are analogues of the three operators D a , K a , and M a we defined to derive vector TAM waves, but they act 
on the Hilbert space of all functions of n; i.e. they do not act on the radial coordinate r. This new set of operators 
satisfies precisely the same algebra as the set {D a , K a , M a }. They are orthogonal to each other, 

N a K a = K a N a = 0, JV a MJ = 0, M ±a N a = 2, K a M^ = M± a K a = 0, (40) 

and they are normalized to 

(N a )^ N a = N a N a = 1, (K a ) f K a = -K a K a = L 2 , (M± a )^ M" = — M]_ MJ = L 2 . (41) 

As operators in the Hilbert space, their hermitian conjugates are 

(N a ) f =N a , {K a )* = -K a , (M ±a ) f = -M ±a + 27V a . (42) 

Furthermore, they satisfy the algebraic relations, 

[N a , N b ] = 0, [Mxa, N b ] = {g ab - N a N b ) , [M± a , M xb \ = -e ab cK c , 

[K a , K b ] = e abc K c , [K a , N b ] = e abc N c , [K a , M ±b ] = e abc M c ± . (43) 

The two operators K a and M± a generate the two transverse- vector spherical harmonics Vjj M j tt (n) and Y(jm)u(**)> 
in terms of E/B modes, while N a generates the longitudinal vector spherical harmonic Yh M s (n) in the normal 
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direction. In summary, 

^(jM)a( ft ) = -j=====K a Y {J M){n) = iY ( J JM)a (n), 



(44) 



normalized so that 



an y, ' n 



(jm) F (J'M> ( fi ) = SjJ'S M M'S, 



»a/J, (45) 



for a,/3 = {-E, £>, L}. Note that although the relations between the OAM and E/B/L vector-sphcrical-harmonic 
bases in Eq. (|44| resemble those between the OAM and E/B/L bases for TAM waves in Eqs. (|35|) and (|37|) . there are 
subtle, and important, sign differences. The minus sign in the definition of N a is chosen to match D a = —k a . The 
sign convention for E/B vector spherical harmonics is chosen so that a rotation of the E mode by +90° about the 
outward normal direction (i.e. the direction of n) yields the B mode. 

Finally, we can write the E/B/L TAM waves in terms of the E/B/L spherical harmonics: 

= Jj(kr)Y ( % I)a (h), 

*(&).(*) = 1 + ]J ^P] F (W fi ) - *VWTT)^Y { L JM)a (n), 

= -^^(^^^aln)-*^!^)^^!!). (46) 

Although the mode functions arc orthonormal, we now see that they are not orthogonal at each point. Although the 
L and B modes are everywhere perpendicular and the E and B modes everywhere perpendicular, the L and E vector 
TAM waves are not always perpendicular. The B mode has components only in the 8-<f> plane; i.e., n a ^>^ JM ^ a (x) = 0. 
The E and L modes most generally have components in the tangential plane and along the normal n a . In the far-field 
limit kr — > oo, however, the three modes are asymptotically perpendicular to each other. 

6. The plane wave expansion for vector fields 



We now determine the transformation between the vector plane-wave basis and the vector TAM-wave bases. We 
start with the OAM basis. Since the ^ , |jA/)a( x ) constitute a complete basis, we may write, 

e (k)e ik '* = ]T ivi l A\ JM) (k)* l * M)a (x) = J2 4™% M) (k)j^r)Y/ JM)a (n). (47) 

I.JM I.JM 

Here i a is a (unit) polarization vector for the wave. The coefficients A^ JM ^{k) may be obtained by writing e lkx in 
the usual scalar plane- wave expansion, 

e (k)e ik ' x = ^47ri i iK^)y ( L)(k)y- (M (fi)£a(k). (48) 

lm 

We then use orthonormality of the Y, l JM ^ a (h) to infer that 



We can similarly expand in terms of L, E, B modes, or helicity modes, as 



4/M)(k) = £ a (k)r// M)a (k). (49) 



e a (k)e ik ~ = J2 E 47ri-/ ^M)(k)*(} a M)o (x), i a (k)e^= E 47riJ ^M)(k)*f} A M)a (x), (50) 



a=L,E,B JM A=-1.0,l JM 
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in terms of expansion coefficients 

\JM) 



A?JW& = ia M Y u*M)a(n 4jM)(k) = i a (k)Y ( Y M)a (k) 



Here, the spin-1 vector spherical harmonics are Y,\^\ = 2 1 I 2 



spin-1 spherical harmonics x^jMl(n) [26j,|27j by 



! (JM)a 

- (JM)a 



yE —j— jyB 
I (JM)a ^ 11 (JM)a 



f)»W- ( 51 ) 

These are related to the usual 



^,(k)YA (k) = _ A Y (JM) (k)5 AV , A, A' = 0,±1. 



(52) 



Here, e^(k), for A = 0, ±1, are the polarization vectors for a vector plane wave with wavevector k and helicity A. This 
equation defines our phase convention for i\. In terms of basis vectors in spherical coordinates, these are defined as 



£ =n , 



1 



■■±i -T^=(f"Ti0 a )- 



(53) 



7. Expansion of vector fields and power spectra 
An arbitrary vector field V a (x) can be expanded in the OAM basis by 



K(X) =£ £ J W 

JM l=J-l, J,J+1 ° y ' 



k 2 dk 



in terms of expansion coefficients, 



V{ JM) (k)= /d 3 x^ Q (x) W^ M)a (x) 



As we show in Ref. [30], these can also be written as vector-spherical-harmonic transforms, 

V l {JM) {k)= f d 2 kV a (k)Yl; M)a (k), 



(54) 



(55) 



(56) 



of the vector Fourier coefficients. Analogous relations hold for the L/E/B and helicity bases as well. 

For the L/E/B basis, we may use the operator approach discussed above to rewrite the expansion coefficients in 
terms of scalar TAM waves by integrating by parts: 



V ( a j M )(k) 



^x^(x) [^ttft^Cx) 



rf 3 x^(x) 47Ti J 0^ JM) (x) 



d 3 x 



(57) 



where O" = {D a , K a , M a }, and the hermitian conjugates of each operator are as given in Eq. (|32|) . Explicit expressions 
for the expansion coefficients are 

r r "I * 

(58) 
(59) 
(60) 





j d 3 x 




^a^Cx), 


V(JM){k) = 


j d 3 x 


47ri J *(> M) (x) 


*(-A' Q )^(x), 


V ( E JM){k) = 


j d 3 x 




*(-M a + 2 J D a )F a (x 



Likewise the coefficients for the helicity basis are 



v ( TUk) = v ( L JM) (k), vggfc*) 



V2 



V { E JM )(k)TiV( 



JM) 



(61) 
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In other words, the expansion coefficients for vector TAM waves are the same as the coefficients of the scalar TAM 
waves for following three scalar functions: 

V L (x) = D a V a (x) = ^V a y a (x) = iv ■ V(x), (62) 
V B (x) = -J( a f(x)= eak i'VT(x) = [xxV].V(4 (63) 
V B (x) = {-M a + 2D a )V a (x) = -p abc V b K c V a (x) + 2^-V a y Q (x) = % - [{V x (x x V)} + 2V] • V(x). (64) 

These scalars may be useful to calculate the theoretical expectation for TAM-wave coefficients. 

Suppose now that V a (x) is written in terms of its longitudinal and transverse parts and that these have power 
spectra -Pi(fc) and Pr(k), as defined in Section TlV Al It follows then that 

(f^W*)]' = PT(k)Sjj'S MW S a p^-5 D (k-k'), (65) 
for {a, (3} = {E, B}. Similarly, 

([VMW]*Vfr M , ) (k')) = PL{k)5jj,5 M M> { ^5 D {k~k% (66) 

for the longitudinal modes. The projections of a linearly-polarized transverse-vector plane wave onto the E and B 
vector TAM waves have equal amplitudes. Therefore, the power spectra for the E and B modes must always be the 
same for a realization of a statistically homogeneous random field. 

V. SYMMETRIC TENSOR FIELDS 

A. Introduction and Plane Waves 

We now consider solutions to the Helmholtz equation, (V 2 + fc 2 ) /i a f>(x) = 0, for a symmetric tensor field h a b{x) = 
h(ab)( x ) = [habfe) + ^6a( x )] /2. The most general such tensor field can be decomposed into a trace component h(x), 
a longitudinal component £(x), two vector components w a (with V a u> a = 0), and two transverse-traceless tensor 
components h^ b (which satisfy V a /i^ b T = and h a a = 0), as 

h ab = hg ab + fv o V 6 - JfahV 2 ^ £ + V (a w fc) + h T J. (67) 

Our goal is to derive rank-2 tensor solutions to the Helmholtz equation, of definite total angular momentum, for these 
different components. 

We begin, though, by reviewing the Fourier decomposition of the rank-2 tensor field. Each Fourier component 
of the tensor field can be expanded as h a b(k) = ^2, s £* b (k)/i s (k) in terms of six polarization states e* 6 (k), where 
s = {0, z, x, y, +, x}, for the trace, longitudinal, two vector, and two transverse-traceless polarizations, respectively, 
with amplitudes h s (k) [28|. The polarization tensors satisfy e sab e a \ b = 2S 8S >. The trace polarization tensor is e° ab oc S a b, 
and the longitudinal is e z ab oc (k a k b — k 2 5 a b/3)k~ 2 . The two vector-mode polarization tensors satisfy e* b oc k( a w b j v 
where w* ,v are two orthogonal (w xa "w y = 0) and transverse (k a w a c,y = 0) vectors. The two transverse-traceless 
polarization states have k a eZ b x = 0. 

The two vector (spin-1) modes x, y can alternatively be written in terms of a helicity basis by defining two hclicity-1 
polarization tensors e^ b = (e* b ± is^ b )/y/2. Similarly, the two transverse-traceless (spin-2) modes +, x can alterna- 
tively be written in terms of a helicity basis by defining two helicity-2 polarization tensors e^ b = (e^ b ± ie* b )/y/2. 

In general relativity, power spectra Ph{k) for gravitational waves (transverse-traceless tensor fields) h a b are defined, 
for example, by 

(A.(k)Mk')) = ^(27r) 3 Mk-k')^^, (68) 

for s, s' = {+, x}, so that 



(h ab (k)h ab (k)) = (2TrfP h (k)6 ss ,S D (k - k'). (69) 
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B. TAM waves 



Our aim now is to find tensor- valued functions h a j,(x), solutions to the tensor Hclmholtz equation for wavenumber 
k, that transform under spatial rotation as representations of order J. These will be cigenfunctions of total angular 
momentum J = L + S. Here the spin can be either S = for the trace of h a \, or S — 2 for the trace-free part. The 
expansion of the trace is simply in terms of scalar TAM waves. We will therefore focus our attention in the following 
on trace-free rank-2 tensors. Therefore, S is now the 5 = 2 spin associated with the vector space spanned by a set of 
basis tensors at each spatial point. 

We start by constructing a rank-2 spherical basis t™ bl for rh = ±2, ±1,0, that transforms under rotations as a 
representation of order 2 by taking direct products of the order-1 spherical basis, 5 

C= J2 (lmilm 2 |2m)erer. (70) 

mim2 

Using orthonormality of Clebsch-Gordan coefficients, these are normalized to (i mi ) ab i™ h 2 * = <5 mim2 . 



1. The orbital- angular-momentum basis 



We begin by expanding the five components of the rank-2 tracclcss tensor in terms of five tensor TAM waves of 
definite orbital-angular-momcntum-squared L 2 for each total angular momentum JM, as 

*(jM)afcW = 31 ( kr ) Y ljM)ab^) = J2(^lm\J M) j, (kr) Y {lm) (n)tf b , I = J - 2, J - 1, J, J + 1, J + 2, (71) 

mm 

an equation that also defines the OAM tensor spherical harmonics Y^ M ^ ab (ii). These OAM tensor spherical harmonics 
of fixed orbital angular momentum satisfy the orthonormality relation, 

J cPhYlf*) (n)F ( Sr MOab (n) = 8 w 8jj,8 MM , . (72) 

The demonstration that the $A ' AI ^ ab (x) constitute a complete basis for traceless symmetric tensors on W 3 , and that 
YrjM)ab constitute a complete basis for three-dimensional traceless tensors on § 2 , are straightforward and similar 
to the analogous proofs for vector harmonics presented in Section HVB II 



2. The longitudinal/vector /transverse-traceless basis 



We now proceed to write the five traceless tensor harmonics for each JM in terms of a longitudinal (L) component, 
two vector components (VE and VB), and two transverse-traceless components (TE and TB). In Appendix IB1 we 
derive the divergence of the tensor spherical waves of fixed orbital angular momentum in terms of vector spherical 
waves to be 



j-i ykj-i 



2J-1 



JM)b 



(x), 



k,J 



.7-1 

2(2.7+1) * (JM)b 

(J+l)(2J+3)~ lT ,fc,J 
6(2J-1)( 2J+1) * {JM)b 

J +2 jvfc.J 
2(2.7+1) *(JM)b 



■7(27-1) 



k,J+l 



6(2,7+l)(2,7+3) (JM)b 



(x), 



(x), 



.7+2 ,T t fc,J+l 
2.7+3 {JM)b 



(x), 



J -2, 
J-l, 
J. 

J+l, 
J + 2. 



(73) 



Note that the divergence of a tensor of fixed total angular momentum JM yields a vector of the same JM, since we 
have acted with V a , an irreducible- vector operator. 



5 Tildcd indices like rh are reserved for the order-2 spherical basis. 
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The transverse-traceless modes. We can, from these results, immediately construct two linear combinations of 
^(jM)a!)( x )' °f different parity, with vanishing divergence, 

^te / (J+l)(J + 2) \ 1/2 j_ 2 / 3(J-l)(J + 2) \ 1 / 2 

V(JM )ab W- { 2 (2J-1)(2J + 1)J W (-/M)a6W 1 V ( 2 J_1)(2J + 3)J V 0"O«&W 

1/2 



/ J(J-l) \ 7 , T ,J +2 

+ \2 (2J + 1) (2 J + 3)J *(JAOaiW' 

*f/M )ab (x) - (^) 1/2 ^^ )ab (x) - (^ r ) 1/2 *(% a& (x) 



(74) 



These two spherical waves form a basis for the transverse-traceless (TT) part of the tensor. We label them E and B 
according to their parity, (— 1) J or (— 1) J+1 , respectively. 

The vector modes. The divergence of the vector component of the rank-2 tensor yields a divergence-free vector field. 
Recalling that the vector harmonics i^^ JM ^ a (x) = ^>^ JM ^ a (x) arc divergence-free, we can construct one vector mode 

of the tensor field by taking the other orthogonal linear combination of ^(jjvf) a f>( x ) an d ^(jM) a &( x ) - Likewise, the 
other vector mode should have a divergence that is proportional to ^fjM) a ( x ^> but it should also be orthogonal to 
the transverse-traceless modes we already obtained. After some algebra, we find the two vector modes of the tensor 
to be 

VE / 2(j-i)(j+i) \ i/2 j_ 2 ( 3 y /2 j 

V(JM)a b W - ( v(2 j_i)(2J + l) ) *(JM)abW + y Jpj — 1) (2 J + 3) ) W 

2j(j+2) y /2 iT „ /+2 



(2J + l)(2J + 3) 



(77) 



These form a basis for the vector part of the tensor field. The divergences of these basis functions are, 

\v a % B M) M = ^ffjM)^ \^ v {J E M)ab {*) = ^ff JM)b {x). (76) 

It then follows that the vector TAM waves of the tensor field can be obtained from the transverse-vector spherical 
waves through, 

*(SM)a b (x) = —£j= (V Q *f JA/)6 (x) + V b *f JA/)a (x)) , *^ )a6 (x) =—A= (Va*fjJl,)6(x) + V b *f JA/)a (x 

The longitudinal mode. The last orthogonal linear combination of the orbital-angular-momentum states, 

ri ( 3(J-1) J V /2 y ,_ 9 / J(J+1) V /2 T J , v 

*( Jtf )'' W " U(2J-1)(2J+1) ) VjOoftW + l v(2 j_i)(2J + 3) J *(-™)a&W 

3 (■/ + !)(■/ + 2) V /2 . rJ+2 (x) (78) 
,2(2J+l)(2J + 3)y ( JM )» tl J ' 1 J 

decries the longitudinal component. To check, we find its divergence to be 

-j;V a *f JA/)ab (x) = *y|*£ /M)fc (x). (79) 
It implies that the longitudinal mode is the only one that has non- vanishing double divergence, 

-Iv a V b *f JA/)a& (x) = y|* (i/M) (x), (80) 
in terms of the scalar spherical wave vE'(jm)- Using this result, it further follows that, 

*fjM)ab(x) = (VaV b - JffabV 2 ") * (JM) (x). (81) 
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3. The longitudinal /vector /transverse-traceless basis in terras of derivative operators 

As seen in Section TlVB 31 the L, B, and E vector waves can be written by applying the vector operators D a , K a , 
and M a , respectively, to scalar TAM waves. Likewise, we have just seen in Eq. (|81|) that the TAM wave for the 
longitudinal component of the tensor field can be written by applying a derivative operator V Q V b — (l/3)<? ab V 2 to 
the scalar spherical wave. We have also seen in Eq. (|TT[) that TAM waves for the vector components of the tensor 
field can be written by taking a symmetrized gradient of the transverse- vector spherical harmonics; i.e., by applying 
the operators V( if b ) and V^Mm, to scalar spherical waves. We now present an operator approach so that we 
have a complete treatment of symmetric traceless tensors in terms of tensor differential operators, including the two 
transverse-traceless tensor modes. 

Using the three vector operators D a , K a , and M a we have proposed as basic building blocks, we construct five 
tensor operators 

T ab = -D a D b + -g ab , T^ b B = D( a K b)l T^ E = D( a M b ), 
Tj B = K [a M b) + M [a K b) + 2D [a K b) , Tj E = M [a M b) - K [a K b) + 2D {a M b) . (82) 

These are irreducible tensors since they are symmetric and traceless, and they commute with V 2 . Therefore, when 
acting on scalar TAM waves ^^ JM y they generate symmetric tensor TAM waves that solve the tensor Hclmholtz 

equation and have the same total angular momentum JM . We have seen that T^ b generates longitudinal mode, and 
T^ b B and T^ b E generate B and E vector modes, respectively. It is straightforward to show that D a Tj b B = D a T^ b E = 

0, so they generate transverse-traceless tensor modes. The parity-odd Tj B generates the B mode ^ tjM} ab {*) (x 

T '<S B *( : ja/)( x )' wmlc tnc parity-even Tj E generates the E mode, *7^ o6 (x) oc T^ b E ^ k {JM) (x). The five operators 
generate five linearly independent tensor modes, because they are orthogonal according to 

(T^Cj T a ' ab = 0, iia^a', for a, a' — L, VB, VE, TB, TE. (83) 
To normalize the tensor spherical waves, we calculate the norms of those five tensor operators to be 

/rpL y rpL,ab _ 2 /jiVBy rpV B ,ab _ frpVE^ rpVE,ab 

_ J(J+1) 

3 2 

(T T S) t T TE,a b = T TB, ab = *(J±W (84) 

To summarize, the decomposition of the traceless symmetric rank-2 tensor into longitudinal, vector, and transverse 
tensor modes is 

-T Q i b * (JM) (x)= ( 2(2J _^(2J + l) j *( 7 /M)a 6 ( X ) 

' j ( 3(J + l)(J + 2) \\ J+2 



J(J + 1) 



(2J- 1) (2J + 3) 



vectorBmode: *|^ /)ab (x) = - J - J + - T^ b B ^ {JM) (x) 



1 2 -t,j~i ( \ , ( j + 2 \ 



n7M)ab(*) + —— * +^ afe (x) 



2j+ij 1 V2j + iy < JM > 



vectorEmode: *f- )a6 (x) = jj^T™ * {JM) (x) = ( ^^1^ ) ' g^fr) 

' 2J(J + 2) \ * j. 
(2J+l)(2J + 3)y < JM ) 




T J ( \ I 2 J ( J + 2) \ - j +2 / \ 
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transverse tensor B mode : *f/M )a6 (x) = - \ / i/ r , o\' , T Ib B ® ( JM) ( x ) 




J — 1 



2J+W 



transverse tensor E mode : ^TjM)ab( x ) = — 



2(J + 2)! Ta6 

3 (J - 1) (J + 2) 
(2J-l)(2J + 3) 



(J + l)(J + 2) 
2 (2J - 1) (2 J + 1) 



J-2 

(JM) 



«„(*) 



(JM) 



ab( X ) 



J(J-l) 



2(2J+l)(2J + 3) 



^J+2 
*(JA/)a6 



X). 

(85) 



The normalizations arc chosen such that 



16tt^ 



<Tx 



(JM) 



^{J'M') ( X ) = &aP&JJ'&MM' 2 p fo(fc " fc'), 



(86) 



where {a, /?} = {L, I^E>, Fi?, T.B, TE 1 }. Again, the five modes are orthogonal as tensor wave functions in the Hilbcrt 
space, but their tensor values at any given point are not necessarily orthogonal, as we will see below. This orthogonality 
does hold asymptotically in the far- field limit kr — > oo. 



4- Summary and helicity basis 

So far we have constructed two sets of TAM-wave bases for symmetric traceless tensors. The OAM basis ^ , (jM)ab( x )' 
where I = J — 2, J— 1, J, J + 1, J + 2 for each JM, are eigenstates of the square of orbital angular momentum L 2 . We 
have also defined a second basis ^ , ™/A/) a b( x ) in terms of a longitudinal mode a = L, two vector modes a = VE, VB, 
and two transverse-traceless tensor modes a = TE, TB, for each JM. We can furthermore construct a helicity basis 
*(jM)a&( x )' denoted by helicity A = ±2, ±1,0, through 

*^M )a6 ( X ) = ^ (*fjM )afe W ±**fj S M)ab (x)) , 
*fjM )a6 ( X ) = ^ (*(SaO«*( X ) ±^(SM)a 6 (x)) , 

*(JM) a6 ( x ) = *(JM) o6 ( X )- (87) 

These are eigenstates of the helicity operator H = S • p, but for the tensor field, (S e ) ab cd = it a ec9bd + Washed, and 
Pa = -(i/k)V a - 

The orbital basis ^ / (jA/)ab( x )' for J = J — 2, J — 1, J, J + 2, J + 2, the longitudinal/ vector /transverse-traceless basis 
* fjM)af,( x )> for Oi = L, VE, VB, TE, TB, and the helicity basis *f JM)o6 (x), for A = 0, +1, -1, +2, -2, are related by 
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unitary transformations, 



U 



M 



3(J-1) J 
2(2,7-l)(2,7+l) 
2(7-l)(.7+l)" 
(2.7-l)(2,7+l) 







V 

/i 





Vo 
/ 



(J+l)(./+2) 
2(2,7-l)(2,7+l) 





l/\/2 
-i/V2 1/V2 



3(7-1)7 
2(2./-l)(2J+l) 
(7-l)(7+l) " 
(2./-l)(2J+l) 

(7-l)(7+l)~ 
(2./-l)(2J+l) 

(J+l)(./+2)~ 
(2,7-l)(2,7+l) 



.7-1 
2,7+1 



J+2 
27+1 






i/V2 



J(J+i) 

(2,7-l)(27+3) 
3 

(2,7-l)(27+3) 



3(J-l)(.7+2)~ 
(2J-l)(2.7+3) 







1/V2 
-i/y/2 1/V2J 





7+2 
27+1 



.7-1 
2,7+1 



(,7+l)(.7+2) 
(2.7-l)(2,7+l) 



.7-1 
2(27+1) 

.7-1 
2(27+1) 

.7+2 
2(2.7+1) 

.7+2 
2(27+1) 



(2.7-l)(2J+3) 
3 

2(2,7-l)(2,7+3) 



2(2,7-l)(27+3) 

3(,7-l)(.7+~ 
2(2,7-l)(2,7+3) 

3(.7-l)(.7+2)~ 
2(2,7-l)(2,7+3) 



3(.7+l)(.7+2) \ 
2 (2.7+1) (27+57 
27(7+2) 
(27+l)(27+3) 







2(27+l)(27+3) 





.7+2 
2(27+1) 

.7+2 
2(2.7+1) 

.7-1 
2(2.7+1) 

.7-1 
2(2.7+1) 



3(7+1) (.7+2) 
2(2,7+1) (2,7+3) 



J(J+2) 
(2,7+l)(2.7+3) 

7(7+2) ~ 
(2,7+l)(2,7+3) 

(2,7+1) (27+3) 

J(J-l) 
(2,7+l)(27+3) / 



5. Projection onto Tensor Spherical Harmonics 

We now describe the projection of three-dimensional traceless tensor TAM waves onto the two-sphere. We will 
begin by reviewing the projection onto the familiar E/B tensor spherical harmonics @ in the 0-(f) space perpendicular 
to n. We will then generalize these two basis tensor spherical harmonics to include three more that will constitute a 
complete orthonormal basis for three-dimensional traceless tensors that live on the two-sphere, parametrized by n. 

The usual E/B tensor spherical harmonics are defined by 



Y iUAs(n) = J J(J+1)(J 2 _ 1)(J + 2) ("V^Vb + \ 9AB V°V C 



Y ( % )AB (n) = ^ 2JiJ+1){ j_ inj + 2) (tB C V c V A + e/V c V B ) Y (JM) (n), (89) 

where here {A, B} = </>}, and is a covariant derivative on the two-sphere, with metric qab = diag(l, sin 2 9) 
and antisymmetric tensor cab- 

The operator approach developed in Section IIV B 51 for vector spherical harmonics can be generalized to tensor 
spherical harmonics. Recall that the operators {iV a , K a , M± a } satisfy the same algebra as the operators {D a , K a , M a } 
do. Therefore, following the same line of reasoning we construct five tensor operators, 

W^ b = -N a N b + ±g ab , Wl B = N [a K b) , W^ b E = N {a M xb) , 
Wj B = K {a M ±b) + M ±{a K b) + 2N ia K b) , W^ E = M Ma M ±b) - K (a K b) + 2N (a M ±b) , (90) 



since they are irreducible tensors. The last two operators are perpendicular to the radial direction, n a Wj b E = 
n a W^ b B = 0, so they must generate the E/B tensor spherical harmonics, n a Wj b E Y^ JM ^ (n) cx ^jf/j^n) and 

^■ a ^Ib B Y{jM){^) ^(7f/)ab(")' respectively. Here ^(jf/) a& (n) and Y^j B ^ ab (h) correspond to the E/B tensor har- 
monics Y^ M ^ AB (n), Y E JM - )AB (h), but under a three-dimensional orthonormal basis, denoted by lower-case indices 
a,b, . . .. In addition, the W^ b E and W^ b B , when acting on ordinary spherical harmonics Y^ JM ^(h), generate VE/VB 
tensor spherical harmonics Y^ E ^ ab (h) and Yjy|^ af) (n), respectively, with components in both the tangential plane 
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and normal. Similarly, W^ b generates the longitudinal tensor spherical harmonics Y^ M ^ Qb (n). For the tensor spherical 
harmonics, the terms "longitudinal," "vector," and "transverse," refer to the nature of the tensor components with 
respect to the normal n a . To be more specific, we define 



^(JM)ab(fi) 



^(JM)ab(fi) 



^W a L b Y (JM) (n), 



Y(jM)ab(**) — 



2(J + 2)! 



(J -2) 



2(J + 2) 



■Wj b B Y iJM) (n). 



(91) 



In terms of the OAM tensor spherical harmonics YK M s ab {n), f° r I = J} JH, J±2, the longitudinal/ vector/transverse 



traceless basis Y^ M ^ ab (h), a = L, VE, VB, TE, TB, are found to be 



longitudinal: Y ( L JM)ab (h) 



vector E mode : Y^ )ab (n) 



+ 



3J(J- 1) 


2(2J-1)(2J 


+ 1) 


J(J + 1) 




(2J-1)(2J + 




2(J-1)(J + 




(2J- 1) (2J + 




3 





(2J- l)(2J + 3)/ (JM) 



3(J + 1) (J + 2) 

2(2J+l)(2J + 3)y J ( JM ) 



2J(J + 2) 



(2J + l)(2J + 3) 



vectorBmodc : 3#& )o6 (fi) = i (j^-M ' Y^Jn) i (±t| j ' Y ( J + M ) ab (n), 



transverse tensor E mode : ^7jM)ob(") = — 



(J + l)(J + 2) 
2(2J- 1) (2J+ 1) 

/ 3 (J - l)(J + 2) 
\(2J- l)(2J + 3) 



J-2 



(JAf)ofcW 
(,7M)ab(fi) ~ 



J(J-l) 



2(2J+l)(2J + 3) 



2 (JM)a6 W' 



transverse tensor B mode : ^jM)ot(°) = 
These are normalized to 



f J + 2 , 



V2J+ 1 



d 2 h 



Y{JM) ( fi ) I F (J'M')a6 ( ft ) = SjJ'^MM'Safi, 



(93) 



for a,/3 = {L,VE,VB,TE,TB}. Note again that although the transformations between the OAM and 
L/VE/VE/TB/TE bases for the tensor spherical harmonics resemble those for the transformations between the 
analogous bases for TAM waves, there are important sign differences. 
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In terms of the tensor spherical harmonics, the TAM waves ^kjM)ab( x ) can ^ e written, 
W = - \ (jj( kr ) + 3 9J (kr)) Y ( L JM)ab (ft) - y/3J(J+l)fj (kr) Y^f I)ab (ft) 



y/3J(J+l)(J-l)(J + 2)jj (kr) vTE 

2 f kr f r (JM)abW> 



( x ) = - V3J(J + l)fj (kr) Y { L JM)ab (ft) - (jj (kr) + 2g., (kr) + 2fj (kr)) Y™ )ab (ft) 
- V(J-l)(J + 2) (/, (fcr) + 2^ (ft) , 



*Jf^ (x) = - V3J(J+1K /' 1)(J + 2 ^^ ( W W - V(J-D(J + 2) [f, (kr) + y- )ab (ft) 

- ^ (-M*r) + 9 J (kr) + ifj (kr) + 6^5?) Y { T } E M)ab (ft) , 
W = " < (iS(fcr) - ^) - ^ -DG/ + 2)^% V (A) , 

(*) = - 'WT^^V (ft) - i (j'Akr) + 2^M) Y™ )ab (ft) . (94) 

Here we have introduced radial profiles, 

fj (x) = ±l^, and 5J (z) = -j^)-2/ J (z) + (J-l)(J + 2)^. (95) 

The components proportional to the transverse Y^jf^ ab (h) and Y^j^ ab (h) harmonics are projections onto the 2- 
sphere. These are the components of principal interest for angular measurements on the sky. Eq. (|M| shows that 
the different tensor TAM waves are not everywhere orthogonal, even though they are orthonormal, although they do 
become asymptotically orthogonal in the kr ^> 1 limit. Both the VB and TB tensor TAM waves have projections 
onto the VB and TB tensor spherical harmonics. The L, VE, and TE TAM waves have projections onto the L, VE, 
and TE tensor spherical harmonics. The phases in our definitions of the tensor spherical harmonics arc chosen so 
that a rotation of a TE (VB) mode by 45° (90°) about n produces a TB (VB) mode. 



6. Plane-wave expansion for traceless tensor fields 



start with the OAM basis. The ^ rjM") a b constitute a complete basis, and we can write 



We now determine the transformation between the tensor plane-wave basis and the tensor TAM-wave bases. We 

l,k 

(JM)c 

l,k 
(J A 

UM 1.7 M 



£ ab (k)e^ = J2 4™ l B{jM)fc)* l (jM)ab(x) = E ^ l B\ JM) (m(kr)Yl JM)ab (n). (96) 



Here e a &(k) is a normalized polarization tensor for the plane wave. We use orthonormality of the Y^j M ^ a (\\) to obtain 
the expansion coefficients, 

B[ JM] (k)^e ab (k)Y { l / M)ab (k). (97) 
We can similarly expand in terms of ^ r'j*M)ab' ^ or a = L, VE, VB, TE, TB modes, or hclicity modes, as 

i ab (k)e^ = ^^4 7 ri^ M) (k)*ff M)o6 (x), e a ,(k)e*- x = £ ^4 7 ri^ M) (k)*^ M)ob (x), (98) 

a JM A=0,±l,±2 JM 

where the expansion coefficients are 

B? JM )(ft = ^(k)^M )Qb (k), 5(W fc ) = iab ( k KM)a b &- (99) 
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Here, the spin-2 tensor spherical harmonics are defined as Y* JA ^ b 



0-1/2 



yTE —J— jvTB 

1 (JM)ab ^ " (JM)ab 



- 1 (JM)ab 



o-l/2 yV£ i ,y7B 
z I 2 (JM)ab 31 ll (JM)ab 

MM by 



and ^/j^^;, = ^fjMlab' These are related to the spin-2 spherical harmonics ;\Y(jM)(n) 



£ ^(k)r ( A /M)Qb (k) = _ A y (JM) (k) fc V) for A, A' = 0, ±1, ±2. 



(100) 



Here, £^ b (k), for A = 0,±1,±2, are the polarization tensors for a tensor plane wave with wavevector k and helicity 
A. This equation defines our phase convention for e^ h . In terms of basis vectors in spherical coordinates, these are 
defined as 



~ab 
£ ±1 



V2 



s ±1 n +e_nn 



±i' 



?ab 
£ ±2 



-£ii£ 



±1 C ±1' 



£n = 



3 (\ 



2 \ 3 



:S a b - h a n 



(101) 



where eg and are defined in Eq. ([53 



7. Expansion of tensor fields and power spectra 
An arbitrary symmetric tracclcss tensor field /i a f,(x) can be expanded in the orbital-angular- momentum basis by 

(102) 



with expansion coefficients, 



JM l = J.J±l,.J±2 J V ' 



h\ JM )(k) = / d 3 x/^(x) [Wtt^^x) 



These can also be written as tensor-spherical-harmonic transforms |3C 



t: 



■l,k 
(JM) 



= I d^T* b (k)Yl; M)ab (k 



(103) 



(104) 



of the tensor Fourier amplitudes. Again, similar relations hold for the L/VE/VB/TE/TB and helicity bases as well. 
The expansion coefficients for the L/VE/VB/TE/TB basis can also be re-written, by integrating by parts, as 



Kj M )(k) = / d 3 x/^(x) 4tt^?}* (x) 



d 3 x^ b (x) A m J T^l JM) {*) 



= / d 3 x 



J.T.fc 



where T" & are the operators defined in Eq. (|82j) . Their hermitian conjugates are given by 



IT. 



L\1 



1 



+ ^9ab, (T^) 1 = -M (a £» 6 ) + 2D a D b: (T^f = -K {a D, 



b), 



(T™^ = M (a M b) - K [a K b) - AM [a D b) - 2D (a M b) + 8D a D b , 
(T^ 5 ) 1 = M (a K b) +K {a M b) -2D {a K b) -AK {a D b) . 



(105) 

(106) 

(107) 
(108) 



We can thus write the expansion coefficients for tensor TAM waves as coefficients of the scalar TAM wave for the 
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following scalar functions: 



-D a D b + -g ab 



h ab (x) 



k 2 



V a V & - -,9abV 2 



1 



^e ac d V c K d V b + ^e bc d V c K d V a + 2V a V t 



h ab (x), 



h v "(x) = --K (a V b) h ab {*), 



(109) 
(110) 
(111) 



h TE (x) 



^\-\ (e a c d V c K d e b JV e K f + e bc d V c K d e a JV e K f ) + 2e ac d V c if d V b + 2e fcc d V c A rf V a 



-V a e hc a V c X d + \7 b e ac a \7 c K d - 8V a V b ^-K (a K b) 

1 i 
~2k 



h ab (x), (112) 

[e ac d V c ^K fc + e bc d V c K d K a + K a e bc d V c K d + K b e ac d V c K d + 2V (a K b) + 4K (a V b) ] /i ab (x). 

(113) 



Here, K a = e abc x b V c . 

If h ab (x) is written in terms of longitudinal/ vector /transverse-traceless parts, and if these have power spectra PlQs), 
Pv(k), and P T (k), then 



([h{jM)(k)] hl J1M1) {k')) = P T (k)6jj,6 M M'6 a0 [ -^-S D (k - k'), a,p = 
( [h?jM)(kj\ * h^, Mt) (k')) - Pv(k)Sjj>6 MM >5 aP ^-6 D {k - k'), a, (3 
( [h(jM)(k)} * hfj, MI) (k')) = P L (k)5jj,S MM ,5 aP ^5 D (k - k'). 



= TE, TB. 



VE, VB. 



(114) 
(115) 
(116) 



For both the vector and transverse-traceless modes, the E modes and B modes have the same power spectra, a 
consequence of statistical homogeneity. 



VI. CALCULATION OF LENSING POWER SPECTRA 



In this Section we provide as an example of the TAM-wave formalism the calculation of lensing power spectra by 
density perturbations and by gravitational waves. We will reproduce results from previous work [3 11 - 4351 ] . which were 
obtained with the Fourier expansion. For clarity, we only take into account the weak-lensing contribution from the 
deflectors (density perturbations or gravitational waves) along the line of sight. However, the measured weak-lensing 
signal also contains other contributions. In particular, metric shear, gravitational-wave effects at the source location 
|3l| , and tidal alignment dominate the power spectrum for lensing by gravitational waves [361 ] . 

Our aim here will be to calculate the lensing deflection field, 



A a (n) = 



n„ 



rjo - V 



drj' 



mi 



n b h eb - - {r}' - ij) h c n d V e h cd 



(117) 



(V,(»7o-»)')fl) 



where n a t, = g ab — n a n b is the projection tensor onto the tangential plane. Thus, the deflection field has no radial 
component and can be viewed as a two-dimensional vector field on the two-sphere. Here, h ab (rj, x) is a (rank-2) tensor 
metric perturbation evaluated at conformal time rj. 

The deflection field A a (n) on the two-sphere can be decomposed into gradient and curl components, 



A„(n) = M ±a <p(h) + K a n(h) 



(118) 



where M± a and K a are the two transverse-vector operators in Eq. (f39|) . and <p(fi) and f2(n) are projected lensing 
potentials. Here we will calculate the contribution A Q (n) that arises from (a) a scalar TAM wave of angular momentum 
JM; (b) an E mode transverse-traceless TAM wave of JM; and (c) a B mode transverse-traceless TAM wave of JM. 
We will then be able to reproduce the power spectra Cj^ and Cj Q for density perturbations and gravitational waves 
that have been obtained earlier by considering individual Fourier modes. These power spectra are defined by, 



( ( P(JM) ( P(J'M')) = SjJ'SmM'C^ 9 , 



JM)&\j>M')) 



1JJ'0MM'Cj Q , 



(119) 
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(120) 



where <P(jm) ancl Q(jm) are spherical-harmonic coefficients for ip(n) and 0(n), respectively. The deflection field can 
be expanded in terms of vector spherical waves by, 

A a (n) = ]T y/J(J + l) [p {JM )Y { E JM)a (ft) + n (JM) Y { B JM)a (n)} . 

JM 

A. Scalar metric perturbation 

Suppose we have a single TAM wave for a density perturbation. This is described by a metric perturbation, 

hM, (m - V')n) = 4^r sca (fc)^M$f/ M) *f JM) (( ?? o - v')n)g ab , (121) 

where Di(rj) is a linear-theory growth factor, a(r/) is the scale factor, and 1S the primordial amplitude of the 

Newtonian potential for wavenumbcr k and total angular momentum JM. We then have 

-\ fa> rf) h b n c v a h bc = -2 (V - „) Ar sca (fc)^^ [v Q *f JM) ((no - v') ft) 

= 2^ (jf - „) At~ (A) ^- k ^f M)a (fo, - ? /) ft) . (122) 

Using Eq. (|46)). which provides the projection of ^,j^(x) onto the plane normal to n a , we find that this single kJM 
mode of the scalar field gives rise to a deflection field, 

A s a ca (n) = y/J(J + l)$f/ M) ^} ca (fc)y ( 5 M)o (n), (123) 

with 

*n*> = fi- <*> £ * (t (w - ,0) ■ (124) 

We thus see that a given TAM wave of JM gives rise only to spherical harmonics in the deflection field of the same 
JM. The absence of a curl (B) mode is a consequence of the fact that the longitudinal TAM wave <tJ' JM j (i (x) has 

no projection onto Y^ M ^(h) [cf. Eq. (|46|) ] . We can equivalently conclude that this particular kJM TAM mode of 

the potential <J> gives rise to a spherical-harmonic coefficient a^ JM ^(k) = ^'f M ^F s j ca (k). The E mode deflection-angle 
power spectrum from the complete random field is then given by summing, 



C 



EE 



J(J+l)Cf^ J k 2 dkP$(k)\Fj(k)\ 2 , (125) 



over all k modes with this JM, in agreement with results obtained by summing over Fourier waves, rather than TAM 
waves. 



B. Tensor Metric Perturbations 

The TAM formalism will have more power, however, for tensor metric perturbations. So consider now a TAM wave, 

h ab (v', (m - V» = h k {J X M) T(k, v')* k ( jM)a b iino - rf)A) , (126) 

of a transverse-traceless metric perturbation. Here, we will take X to be either E or B (although we could have 
alternatively considered A = ±2 modes in the hclicity basis), T(k, rf) gives the time evolution of modes of wavenumbcr 
k, and fyJ^Q is the primordial amplitude of the mode. From Eq. (|117p . we will need to calculate the tangential 
projections of ™ a ^(jM)ab( x ) and n a n c Vf, v I'J C j^ f)ac (x). Using the transformation in Eq. ([88]) between the OAM and 
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L/VE/VB/TE/TB bases and the relations, Eq. (|94|) . between the tensor TAM waves and the tensor spherical 
harmonics, one can show that 



a, j, i 



&3J-2{kr)Y^ )b (n), 
-jj-x(kr)Y^ JM)i (n), 



2(2 J+l) •■ 



(JM)ab 



(x) = 



J(2J-1) y.J+l 



A 7 ( kr ) (_„ I (./+1X2./+3)" J-l (h) 

JW) y Y 6(2.7-l)(2.7+l) l (JM)b\ ll > ^ Y 6(2J+l)(2J+3) 1 (JM)b 
- \/wi7TTj3W ( kr ) Y (JM)b . 



(ft)), 



= J-2, 
= J-l, 

= J 

= J + l, 
= J + 2. 



(127) 



It then follows that, 



™°*(JM)a 6 (*) 



(J 2(2 jj J l| 2 2) (^-i(fcr)+i J+1 (fcr))F ( S M)6 (fi) 



(J-l) (J + 2)jj(fcr) vB 

2 to^ Y (JM)bW, 



(128) 



which lies entirely in the plane normal to fi, and 



(J - 1)(J+1)(J+2 ^i,-2(^ + i^))y ( - (fi) 



2 (2J- 1)^(2 J+l) 



' J(J-l)(J + 2) 
2(2J + 3) 2 (2J+ 1) 



(jj(kr)+ ]J+2 (kr))Y { J ^ )h (h) 



l (J-l)(J + l)(J + 2)jj- 1 (kr) j_ t 



2(2J+1) 



J,J - 1)|J + !),W, '^.» P») 



2(2 J + 1) 



kr 



This vector has components both in the normal direction and in the tangent space. The projection onto the tangent 
space is 



n fcc ™ Q vff JM)ac (x) = - 

The first term in Eq. (|117|) contributes 



(J-l) (J + 2) 1 
2 kr 



j'j (kr) + 



jj (kr) 
kr 



Y( E jM)b(n)- 



(130) 



A GW,X W = / (J mi + V Y( X M)a (ft) / r ' */ 



(770 - r?) fc (770 - 7?') 



7t/j (k(vo ~ V')) , (131) 



where the radial functions are 



fj (kr) 



I ijj (kr) , 
[j'j(kr 



3 J ( kr ) 
kr '■ 



for X = B, 
for X = E. 



(132) 



We now turn to the second term in Eq. (j!17[) . that proportional to ficfi^ e h cd . Using V a 7ib = (g a f, — h a fib) jr = 
n a b/r, we can write 



h c h d V e h cd = V e (n c n d h cd ) - (V e 7i c ) h d h cd - n c (V e h d ) h cd 
= V e (n c n d h cd ) - 2 (V e n c ) n d /i cd 
2 



V e (h c h d h cd ) 



(133) 



(Vo - ?/) 

Given that n ae n e c = H ac , we see that the second term here is similar to what we calculated before. It thus contributes, 

r \, (v'-v)T(k, V f ) 



(J-l)(J + 2) vX 



Y(j M )a («) / dr, 



(r) - 77) k (770 - rfY 



ff(k( m - v ')), 



(134) 
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to the deflection field. Now consider the first term in Eq. (|133|) . We have already seen for the B mode that 
n a ^^ JM ^ ab (x.) is perpendicular to h b . Thus, n a n & \f , ^ M ^ ab (x) = 0, and this term does not contribute to the curl 

(B mode). Using Eq. (fT29f and 



I = J —1, 
I = J, 



(135) 



■y/$frY( JM) (&), l = J+l, 



wc find 



n a n b Vf JM)ab (x) 



(J + 2)! jjjkr) 

— 2-YjM{n). 



2(J-2)! (k r y 



(136) 



When we now take the gradient —M± a = rIi a bV b of this in the tangential plane, the operator acts only on the 
spherical harmonic. Using M±_oY{jm) = \J J(J + ^}YM M s a , we find, 



n ab V»h c h d *Jf M)cd (x) = - 



We thus find that this term contributes 



(J-l)(J + 2) J(J + 1) 



^Jj(k(vo-v'))Y { E JM)a (n). 



A aw, E{3) (fl) = y ( J 1M J + 2 ) j ( j + 1) r ( f M)a (fi) / *, 



/c 2 (770 - 77') 

' , (rf- V )T(k,rf) 



(170 - 77) & 2 (»7o - f/)' 



2 V 2 

In summary, a single B- or E-modc TAM wave contributes a deflection field, 

A GW (fi) = ^ ) yj(jTl)F J GW ' x (A ; )Y- ( 5 Af)a (fi) J 

for X = {E, B}, with 



F GW ' B (fc) 



(J-l)(J + 2) 1 f , T(M0 . , 
2J(J+1) Vo-VJno fe (770 — 77O 



and 



? GW,E, 



(J-l)(J + 2) f^, T(M') 



x{[#(*fa>-V)) + 
The power spectra are then obtained from 



2 J( J + 1) J m k {no - rff 
j.i (k {vo - v')) 



C 



xx 



Hvo - v') 



J{J -i] ].j{k{^-7f)) 

2 770 - 77 fc(?7o - 77') 



F? w ' x (k) 



(137) 



-jj(k(vo-v'))- (138) 



(139) 



(140) 



(141) 



(142) 



for X = E,B, by summing over all k modes with this JM. 



VII. BOLTZMANN EQUATIONS FOR CMB FLUCTUATIONS 

TAM waves can also be used to provide an alternative set of Boltzmann equations to calculate CMB power spectra. 
Our discussion here is preliminary; we leave the full calculation to future work 1181 . There is some overlap, although 
not complete, with what we discuss here and work in Ref. [l7|, and also in Ref. |37l|. 

The radiation perturbation 0(x, qjr/) is most generally a function of position x, the photon direction q, and 
conformal time 77. This perturbation satisfies a Boltzmann equation, a partial differential equation in time, space, 
and in photon direction q. In the standard treatment [IH, one considers a single Fourier mode $(x, 77) = $ke lk x of 
wavevector k of the gravitational potential (or of the gravitational- wave field). The spatial dependence of 0(x, q;?j) 
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must also then be oc e . The q dependence is, however, then expanded in spherical harmonics. Since the end 
result, the power spectrum Cj, is a rotational invariant, one generally then chooses k |j z so that the spherical- 
harmonic expansion for the q dependence of 8k(q;?7) becomes in practice an expansion in Legendre polynomials 
P(cos6 p ) cx Y(./o)(q). 

Alternatively, though, the gravitational potential can be expanded 4>(x, 77) = YlkJM ^(jm)(^)^(./m)( x ) m terms 01 
scalar TAM waves \l/^ JM ^(x) = j.j(kx)Y^j M - j (x) (or for tensor perturbations, in terms of tensor TAM waves). The 
most general radiation perturbation associated with this scalar perturbation can then be expanded in terms of states 
of TAM JM, 

6(x ) q;r 7 )= ]T e^%)s£/ A '(x, q), (143) 

k.JM.U' 

where the total-angular-momentum eigenfunctions (which are also eigenf unctions, of quantum numbers I and I', of x 
and q angular momentum, respectively) are 

S|y/ M (x,q) = ]T (lml'm'\JM}j l (kx)Y (lm) (±)Y (l , m , ) (^). (144) 

mm' 

It now follows that the angular dependence of the observed radiation from a spherical wave with quantum numbers 
kJM will be proportional to Y^ JM ^(q). We take the observer to be at the origin. We then note that the radial 
eigenfunctions ji(kr) all vanish at the origin unless I = 0. Thus, 

6(x = 0, q; V ) = £ ef/^/™(x = 0, q) = 9 Y M (r/)3^ M (x = 0, q) 
w 

= e k 'f M (OOJAf I JM) F ( oo)(x)io(0)y (J M)(q)- (145) 

In the TAM approach, therefore, calculation of the CMB temperature fluctuation boils down to calculation of 
Q Q 'j M (i])- The Boltzmann equation for this particular coefficient, however, will be coupled to those for all 0f/, JM . 
We thus trade the infinite tower of equations for the I = 0, 1, 2, . . . coefficients @i(k) for each wavenumber k for an 
infinite tower I, V — 0, 1, 2, . . . for the coefficients Q,,', J for a particular J. The advantage, though, is that each TAM 
wave of JM contributes only to Cj. Thus, the power spectrum Cj can be evaluated for a single J. There may also 
be conceptual advantages to this approach, even if there are no immediate numerical advantages. 



VIII. CONCLUSIONS 



In this paper we have obtained complete sets of basis functions, specified by their total angular momentum JM, 
for scalar, vector, and tensor fields on R 3 . We have written three such sets of basis functions, one in terms of orbital- 
angular-momentum states, one in terms of an L/E/B or L/VE/VB/TE/TB decomposition of the vector and tensor 
fields, and a third in terms of states of definite helicity. In the process, we have also shown how all five components 
of a rank-2 traceless symmetric tensor field, including the transverse-traceless components, can be written in terms 
of derivative operators acting on scalar fields, a result that may be useful for basis functions beyond those, based on 
spherical coordinates, that we have derived here. 

We have shown how the projections of these three-dimensional vector and tensor fields onto the two-sphere yield 
the familiar E/B vector and tensor spherical harmonics. We found that an E mode on the two-sphere may arise from 
either a longitudinal vector or tensor mode or from E mode vector or tensor TAM waves. Conversely a B mode on 
the two-sphere is seen to arise from a projection of a vector or tensor B mode. We also generalized the two usual 
E/B tensor spherical harmonics to account for the three other possible polarizations of a traceless three-dimensional 
tensor field. We showed how the five TAM waves project onto these five tensor spherical harmonics. 

A realization of a random scalar, vector, or tensor field is usually described as a collection of plane waves with 
amplitudes selected from some distribution. We have shown, however, that a random field can also be realized as a 
collection of TAM waves, and we have shown how the power spectra for these TAM- wave amplitudes are related to 
the power spectra for the more familiar plane-wave amplitudes. The advantage of TAM waves over the simpler but 
more naive outer product of the tensor spherical harmonics with radial wave functions is that such basis functions 
are not necessarily eigenfunctions of the Laplacian. They therefore do not follow a simple evolution during the linear 
regime at late times, and they are not normal modes during inflation. 

The utility of TAM waves in cosmology is apparent given that most observations are performed on a spherical 
sky. Many calculations of cosmological observables, which are usually performed by considering the projection of a 
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single Fourier mode onto a spherical sky, can be performed alternatively by considering a single TAM wave. The 
angular dependence of any observable on a TAM wave of total angular momentum JM must then be a scalar, 
vector, or tensor harmonic (depending on the observable) of that same JM. We showed, as one example, how the 
calculation of power spectra for the lensing-deflection field for gravitational lensing by density perturbations and 
gravitational waves is carried out in the TAM formalism, and we made preliminary remarks on the possible utility of 
the TAM formalism in numerical evaluation of CMB power spectra. The full power of the TAM formalism will be 
manifest most clearly, though, in the calculation of higher-order correlations (e.g., angular bispectra) in models with 
non-Gaussianity, particularly those involving vector and/or tensor fields. The basic idea here is that the Wigncr- 
Eckart theorem guarantees that angular correlations of three TAM waves must be proportional to a Clebsch-Gordan 
coefficient, along with some prefactor that will depend on the tensorial nature of the waves. This will be presented in 
Ref. [IJ. 

The development of the TAM-wave formalism requires considerable technical detail. However, once completed, 
understood, and mastered, it may facilitate the calculation of many cosmological observables. 
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Appendix A: Divergence of the vector harmonics 
1. Gradient of scalar TAM waves 



First, we calculate the gradient, 

Va*f JA/) (x) = V a]J {kr)Y (JM) (h), (Al) 

of scalar TAM waves. Our derivation is based on the Fourier transform. From Rayleigh's formula [Eq. ||5J], we find 
the Fourier transform of scalar harmonics jj(kr)YrjM) (n) as 



^(jM)(q), 



and re- write the gradient as an inverse Fourier transform: 

v a *f, 7M) (x) = Vo y ^ 



9 2/ , v / ^(g~fc) 

2ir [-1) ^ Y (JM)W 



£ 



i-j+i 



q 2 dq^ji(qr)5 D (q - k) 



d 2 qq a Y {JM) (q)Y* lm) (q) 



The radial integral trivially reads 

J 1 2d 9^2.ii(q r ) S D(q -k) = kji(kr), 
and the angular integral can be written as a Gaunt integral by decomposing 

£y (m) (n)(-lfer\ 



Y { i m )(n). 



(A2) 



(A3) 



(A4) 



(A5) 



Then, 



J rf 2 qq a y (JM) (q)y ( t m) (q) = ^-^(-l)^-^ J d'qY^ (q)F (JM) (q)F ( | TO) (q) 



I An 



m+m — m 



3(2J + 1)(2J + 1) / 1 J I 



-in 







£(-i) A V 



3(2J+1)(2/ + 1) / 1 J I 
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1 J I 
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(A6) 
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Here, 



'1 h I3 



1 



= (-1) 



h — l2+m 3 



(lim 1 l 2 m2\hrnz) 



m 1 to 2 -m 3 J ^/2l 3 + 1 
is the Wigner-3j symbol. Combining all, we find the gradient of the scalar TAM wave to be 



(A7) 



! Ill 



*Z) i '~ J+1 J«(*0^|^ (*010|J0) E>a™ (1 - *nZm|JM) K (lm) (fi) 



where we have used that the relevant Clebsch-Gordan coefficient is non-zero only for 

<Z010|J0) = 



i±L 7=7-1 

2Z+1' ' 



j 7=7-1-1 



(A8) 



(A9) 



2. Divergence of vector TAM waves 

We can now derive the divergence, given in Eq. (|27|) . of the vector TAM waves. We start with 



jt(kr) J2 (lmlm\JM)Y {lm) (n)ef 



= ]T (imm\JM)re™V a 3i{kr)Y {lm) {iL) 

m,fh 

and then from Eq. (|A8p . we have 

^ a *;; 7 fc M)Q (x) = fcr^(lmH«>< 
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-WT\ (x) + W— — (x) 



(A10) 



(All) 



The sum can be simplified as 
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m to — Af / \ to to —to' 



= -Jt±i(fer) E ^^J+Tv^m ^ 1 ^"'' Y (l±Uml) (&) = -ji±i(kr)^j ^^Sj, l±1 Y (l±1M) (n), (A12) 
from which follows Eq. (|27[) . 
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Appendix B: Divergence of tensor harmonics 



In this Appendix we return to the use of our usual index notation for vectors and tensors so that vectors and tensors 
can be distinguished by the number of indices. The divergence of the tensor TAM waves is, 

lv a *{jM)ab(x) = J2(2mlm\JM)± {V a 3l (kr)Y {lm) (x)) i* = J2(2mlm\ JM>*f,-)(x)^. (Bl) 

mm mm 

From Eq. (|A8p , we calculate 

*fc)C= E E (l^W|Zm)% m0 (x)<lm 1 lm a |2m)e*«(e* 1 e*») 

= E E C" 1 )*^. lm 2 |2m){lmrm'|Zm}* (iW) (x)er. (B2) 

m2 m'm 

We first work out the sums over to, to, and to. We trade Clcbsch-Gordan coefficients for Wigncr-3j symbols, 
(-l)™(2mlm\JM)(l,-m, lm 2 \2m) {Irhl 1 m' \lm) 
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where we have used the definition, 
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for the Wigner-6j symbol. We are then left with 
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We evaluate the 6j-symbols explicitly and find 
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Note that we automatically obtain eigenfunctions of total angular momentum as a consequence of acting with V a , an 
irreducible-vector operator, on a tensor eigenfunction of total angular momentum. 
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Appendix C: Irreducible tensors 

In this Appendix, we show that irreducible-tensor operators, when acting on a wavefunction, conserve the total 
angular momentum, even though the spin might change. To be more precise, let us consider the spherical harmonics 
Y(jjv/)(fi)i which are eigenfunctions of orbital angular momentum L 2 and L z , for given JM; i.e., 

L 2 F ( JM) = J ( J + 1) Y {JM) , L Z Y (JM) = MY (JM) . (CI) 

Assume we have a group of irreducible-tensor operators O m , for m = + 1,1, that transform as a 

representation of order / under rotations. There are 21 + 1 such operators, and O m YrjM) is a spin-/ object, a tensor 
wavefunction of higher rank. There will be spin operators S a , for a = x,y, z, that act on such spin-/ objects. The 
total angular momentum J a = L a + S a is then defined as the sum of the orbital angular momentum and the spin. We 
would like to prove that 

3 2 O l m Y iJM) = J ( J + 1) G l m Y(j M ) , JAYm = MO l m Y (JM) . (C2) 

Consider a rotation TZ = e 40 Ja acting on the Hilbert space of spin-/ wavefunctions, where O a parametrize rotation 
angles. The orbital angular momentum L a generates rotations of configuration space, and the spin S a generates 
rotations of the internal tensor space; i.e. mixing of the tensor components. Using the fact that L a and S a commute, 
we have 

je a j a fni v J® a L a ( A@ b s b/ ni \ a -ie c L c je d L dV fr ,n\ 

e U m Y (JM)=e "U m je "e d Y( JM y (C3j 

First we note that spin operators rotate O l m contravariantly by acting on the tensor index m (if the tensor basis 
transforms covariantly) , 

e ie°s a0 i n = J20 l mi V l * mi (K), (C4) 

mi 

where we have introduced the Wigner D rotation matrices T> m m {TV)- Meanwhile, the orbital angular momentum 
rotates K(jm) covariantly, 

e m * L «Y {JM) = J2 Y [JMl) V J MlM (K) . (C5) 

Ah 

Finally, O l m , being irreducible-tensor operators, are rotated by the orbital-angular- momentum operators according to, 

e^ L "O l mi e-^ bL > = O l m2 V l m m ■ (C6) 

II) 2 

Prom the above three equations, plus the unitarity of Wigner P-matrices, 

J2 V l n \ mi (TZ) V l m2mi (K) - 6 mm2 , (C7) 

mi 

we find 

e ie°J« i mY{JU) = O l m Y {JMl) V{ hM (K) . (C8) 

Mi 

The derivation holds for any rotation TZ, so we conclude that spin-/ objects O l m Y(jM) transform as a representation 
of order J under rotations, and hence must be eigenfunctions of total angular momentum, as described in Eq. (|C2[) . 

The proof is easily generalized if the Y(jm) arc replaced by spherical harmonics of higher spin. Then the orbital 
angular momentum L a is replaced by total angular momentum, while J a will be the new total angular momentum 
which is obtained by adding the additional spin carried by irreducible-tensor operators O m . 
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